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Abstract 

In the recent paper by one of the authors (MBS) and A. A. Malykh 
on the classification of second-order PDEs with four independent vari- 
ables that possess partner symmetries [TJ, mixed heavenly equation 
and Husain equation appear as closely related canonical equations 
admitting partner symmetries. Here for the mixed heavenly equa- 
tion and Husain equation, formulated in a two-component form, we 
present recursion operators, Lax pairs of Olver-Ibragimov-Shabat type 
and discover their Lagrangians, symplectic and bi-Hamiltonian struc- 
ture. We obtain all point and second-order symmetries, integrals and 
bi-Hamiltonian representations of these systems and their symmetry 
flows together with infinite hierarchies of nonlocal higher symmetries. 



1 Introduction 

In the recent paper [TJ, one of the authors (MBS) and A. A. Malykh ob- 
tained, up to a change of notation for independent variables, the general form 
of second-order partial differential equations (PDEs) with four independent 
variables t, x, y, z, that possess partner symmetries [2H5] and contain only 
second derivatives of the unknown u: 

F = a\{u t yU xz — u tz u xy ) + a2{u tx Ut y — u u u X y) + a^(u ty u xx — u tx u xy ) 
+ a 4 (ut x u tz - u tt u xz ) + a 5 (u tz u xx - u tx u xz ) + a 6 (u t tu xx - u 2 tx ) (1.1) 
+ biu xy + b 2 u ty + b 3 u xz + b A u tz + b 5 u tt + 2b 6 u tx + b 7 u xx + b = 0, 
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with constant coefficients Oj and 6j. Partner symmetries, that make it possible 
to obtain noninvariant solutions of PDEs of the form (II. ip . are generated by 
the recursion relation: 

<Pt = - (a2U ty + a 4 u tz - a 6 u tx + h~ oJoj^t - 

(a 3 u ty + a 5 u tz + a 6 u tt + b 7 )(p x + (axUtz + a 2 u tt + a 3 u tx - bx)(p y 

+ ( - a x u t y + a 4 u tt + a 5 u tx - b 3 )<p z , 
<p x = -{a%u xy + a 4 u xz - a G u xx - b 5 )(p t - (1.2) 
(a 3 u xy + a 5 u xz + a 6 u tx - b 6 - u )ip x + (aiu xz + a 2 u tx + a 3 u xx + b 2 )^ y 

+ ( - a x u xy + a 4 u tx + a 5 u xx + h)ip z , 

where ip and (p are symmetry characteristics [B] and Uo is a constant. In 
( II. ip and (II. 2p . subscripts denote partial derivatives. The transformation 
( II. 2p maps any symmetry (p of the equation (II .ip into its partner symmetry 
(p. 

In jl], we also listed canonical forms to which the general form ( II .ip can 
be reduced by point and Legendre transformations. Among these forms we 
find, along with the first and second heavenly equations of Plebanski [7], a 
new equation that looks, up to a point, like the combination of these two 
equations, which we called mixed heavenly equation: 

Ut y u xz — u tz u xy + u tt u xx — u\ x = e, (1-3) 

where e = ±1. Recursion relation (II. 2p for symmetries of the equation (II. 3p 
becomes 

(f>t = (U tx + UJ())<Pt - Utt(p x + U tz tPy - U t yip zj 

<p x = u xx <p t - (u tx - u )<p x + u xz ip y - u xy ip z . (1.4) 

Note that in our classification heavenly equations of Plebanski belong to 
equivalence classes different from the one to which the mixed heavenly equa- 
tion belongs, that is, they cannot be related neither by point nor by Legendre 
transformations^ 

Another form of a canonical equation from the same equivalence class 
coincides, at e — +1, with the Husain heavenly equation: 

u ty u xz - u tz u xy + u tt + eu xx = 0, (1.5) 

t Quite recently a different classification of integrable PDEs of Plebanski type was done 
in paper [8]. 
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which is an alternative form of a basic self-dual gravity equation arising in 
the chiral model approach to self-dual gravity [5|[TU]. Recursion relation (}1.2j) 
for symmetries of equation (11.51) takes the form 



Though equation (11.51) can be obtained from the mixed heavenly equation 
(11.31) by Legendre transformation ([7.2]) . the main objects of the Hamilto- 
nian formulation of Husain equation, like Lagrangian, symplectic two-form, 
Hamiltonian operators and Hamiltonian densities cannot be obtained that 
way. Therefore, we study Lax representation, symplectic and Hamiltonian 
structures of Eq. (jl.5p independently of those of Eq. ( II. 3p . 

In this paper, we consider mixed heavenly equation and Husain equa- 
tion in a two-component form, which enables us to rewrite the correspond- 
ing recursion relation as a single 2x2 matrix-differential equation and in- 
troduce naturally recursion operators R. Together with the operator A of 
the symmetry condition, which determines symmetries in a two-component 
form, the two operators R and A form Lax pair of Olver-Ibragimov-Shabat 
type [LTIEEI • We construct symplectic and Hamiltonian operators and corre- 
sponding Hamiltonian densities for the two-component mixed heavenly sys- 
tem and Husain system. The proof of Jacobi identity is a simple check that 
the corresponding symplectic two-forms are closed. Thus, both systems are 
set into Hamiltonian form. Applying the recursion operator R to first Hamil- 
tonian operator, for each system, we generate explicitly second Hamiltonian 
operators for both systems and, thus, we show that they are bi-Hamiltonian 
systems, which are integrable in the sense of Magri [13] . 

In section [21 we derive recursion operator for symmetries and Lax repre- 
sentation for mixed heavenly system. 

In section El we present Lagrangian, symplectic two-form, symplectic op- 
erator and Hamiltonian representation of the mixed heavenly system. 

In section HJ using the recursion operator, we obtain explicitly second 
and third Hamiltonian representations of our system and prove that the first 
two Hamiltonian operators are compatible, i.e. they form Poisson pencil. 
Thus, we show that mixed heavenly equation in a two-component form is an 
integrable bi-Hamiltonian system in the sense of Magri. 

In section [5[ we present all point symmetries and second-order symme- 
tries of the mixed heavenly system and calculate a table of commutators of 



(fit = U tz (f y - Uty(fz - £<fx + W V?i, 

fix = Uxz^y - U xy (p z + <f t + U ip x . 
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symmetry generators. In subsection 15.11 we study the action of the recur- 
sion operator on these symmetries. In subsection 15. 2[ using inverse Noether 
theorem for variational symmetries, we determine conserved densities corre- 
sponding to them, which serve as Hamiltonians of the symmetry flows. In 
subsection I5.3[ we study recursions of Hamiltonians for symmetry flows of 
the mixed heavenly system. 

In section we study hierarchies of mixed heavenly system and its sym- 
metry flows. In contrast to our previous studies of bi-Hamiltonian struc- 
tures of Plebanski second heavenly equation [T3] and complex Monge- Ampere 
equation [15] , the mixed heavenly system, though being bi-Hamiltonian, does 
not possess an infinite hierarchy but its hierarchy consists only of two mem- 
bers (the same is true for Husain system). However, some of its variational 
symmetries do form an infinite hierarchy of commuting flows, that contains 
higher nonlocal flows. In particular, we present explicitly the first Hamil- 
tonian flow generating a nonlocal symmetry of our system. We also obtain 
bi-Hamiltonian representations for all local variational symmetry flows. 

In section [3, we convert Husain equation in a two-component form and 
present a Lagrangian, suitable for deriving Hamiltonian form of the two- 
component Husain system. 

In section [HJ we construct a symplectic two-form, symplectic and Hamilto- 
nian operators and Hamiltonian density and, thus, obtain Hamiltonian form 
of the Husain system. 

In section [21 we discover recursion operator and obtain Lax representation 
for Husain system. 

In section [TO], we compute second Hamiltonian operator and prove that 
the two Hamiltonian operators form Poisson pencil. Thus, we conclude that 
Husain equation in a two-component form is an integrable bi-Hamiltonian 
system in the sense of Magri [T3] . 

In section [11], we present all point and second-order symmetries of Hu- 
sain system and calculate a table of commutators of symmetry generators. 
In subsection 111.11 we study the action of recursion operator on symmetries 
of Husain system. In subsection 111.21 we find Hamiltonian densities of varia- 
tional symmetry flows which are conserved densities of the Husain system. In 
subsection |11.3[ we study the action of recursion operator on Hamiltonians 
of the symmetry flows. 

In section [121 we show that the hierarchy of Husain system consists of two 
members only. We obtain an infinite hierarchy of commuting symmetry flows 
that contains higher nonlocal flows. We present explicitly the first nonlocal 
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Hamiltonian flow in this hierarchy. We obtain bi-Hamiltonian representations 
for all local variational symmetry flows of Husain system. 



2 Recursion operator and Lax representation 
of mixed heavenly system 

By choosing ut = v as the second unknown, we present mixed heavenly 
equation (11. 3p in the form of a two-component evolution system 

u t = v, v t = —(vl + v z u xy -v y u xz + e)=Q, (2.1) 

11>xx 

which we shall call mixed heavenly system. Lie groups of symmetry transfor- 
mations of system (12 .ip in the canonical form, when only dependent variables 
are transformed, are determined by the Lie equations 

u T = ip, v T = ip, (2.2) 

where r is the group parameter. The symmetry condition amounts to com- 
patibility of Lie equations (12. 2p and equations (I2.ip : u tT — u rt = and 
v tr — v Tt = 0. We introduce a two-component symmetry characteristic of sys- 
tem (12 $ = ( ^ ^. Then the symmetry condition results in the linear 
matrix equation 

A($) = 0, (2.3) 
where A is the Frechet derivative of the flow (12. ip 



.4 



D± -1 



Q n 2 v z n n i v v n n n x n i u xz n Ux v n 

11 x 11 J - y x J - y y T q, J-^x^zi - L -'t 7/ ^x T 7/ J-'y 7/ 



(2.4) 

Here -D*, D^, D z are operators of total derivatives with respect to t, x, 
y, z. In particular, the first row of (12. 3p yields (ft = ip. Using this relation 
and a similar one for the partner symmetry, ip — (ft, we rewrite the recursion 
relation ( 11. 4ft with wo = and Ut = v, Vt = Q in the two-component form 

ip = ( - QD X + v z D y - v y D z )ip + v x ip, 
<Px= {- v x D x + u xz D y - u xy D z )ip + u xx ip. (2.5) 
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After integrating the second equation (12.51) with respect to x at constant y, z 

( <p \ 

and t and using the notation $ = r I , the recursion relation takes the 

\V> / 

matrix form $ = -R($), where the recursion operator R is the 2x2 matrix 

r = ( VxDx + UxzE> y ~~ u *y Dz ) D x lu ** \ ^ 2 g x 

V -QD X + v z D y - v y D z v x J 

For the commutator of the recursion operator R and operator A of the 
symmetry condition ( 12 ,3j) . computed without using the equations of motion, 
we obtain 

M]= (2-7) 

(v t -Q) x - D- 1 [(v t -Q) xx _n-V„ _^ 

+(«* - «)„£> tf - («« - v) xy D z ] UxKUt v)xx 



v )xy 



[-Q(u t - v) xx + v z (u t 
-v y {u t - v) xz + 2v x {v t - Q) x - u xz {v t - Q) y ~( v t ~ Q)x 
+u xy (v t - Q) Z ]D X - (v t - Q) z D y + (v t - Q)D Z 

Therefore, on solutions of the system (12. ID operators R and A commute and 
hence R acting on any two-component symmetry generates again a symmetry. 
This proves that R is indeed a recursion operator. Moreover, vanishing of 
the commutator (12. 7p reproduces (12.11) and hence operators R and A form 
a Lax pair of the Olver-Ibragimov-Shabat type [THEE] f° r mrx ed heavenly 
system (12. ip . 

3 Lagrangian, symplectic and Hamiltonian 
structure of mixed heavenly system 

We start with the Lagrangian for the mixed heavenly system (12 .ip 

vu t ~\ v ) Uxx + \ u ti. u v u xz ~ u z u xy ) + eu, (3.1) 
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which yields the canonical momenta 

dL 1 . dL , . 

7T„ = — = UWaKc + - [UyU xz - U z U xy ), 7l v = — = 0, (3.2) 

ou t 3 cw f 

that cannot be inverted for the velocities u t and v t , and therefore the La- 
grangian (13. Xp is degenerate. Following Dirac's theory of constraints [16J, we 
treat the definitions (I3.2p as constraints of the second class 

4>u = n u - vu xx - ^ (u y u xz - u z u xy ) = 0, (j) V = n v = 0, (3.3) 
compute the Poisson brackets of the constraints (for details of this procedure 



see 

[<f> i (x,y,z),<p j (x',y',z , )] = Ky, i,j = 1,2, u x = u, u 2 = v (3.4) 
as entries of the 2x2 matrix 

f D x v x -\- v x D x -\- 2 (D z ii X y -\- u X yD z ^ g (Dyii xz -\- u xz Dy) u xx 

x V u xx 

(3.5) 

which is an explicitly skew-symmetric symplectic operator. The correspond- 
ing symplectic two-form is a volume integral Q — J udxdydz of the density 

v 

oj = ^du l A Kij du j (3.6) 
= v x du A du x — u xx du Adv + — [u xy du A du z — u xz du A du y ) . 
The form Q is closed since the exterior differential of ( 13.61) is a total divergence 

du — - du x A du y A du z 

= - (lD x (du A du y A du z ) + D y (du x A du A du z ) + D z (du x A du y A du)\ 
6 (3-7) 

which is equivalent to zero under the volume integral in Q, at appropriate 
boundary conditions. Therefore, Q is indeed a symplectic form and so K, 
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defined by (13. 5p . is indeed a symplectic operator. Hence its inverse is a 
Hamiltonian operator 



Jo = K- 1 



( J_ \ 

1 _l_ 1 y) Vx — ,V' XZ D — ■ 

\ t^xx x Uxx Hxx x l^xx 'Zu X x ^ ^ 



XX 



U xx _ 1 n Uxz_ , U xy n _J_ , 1 n ^xy_ I 



(3.8) 

since it is explicitly skew-symmetric and Jacobi identity is satisfied as a 
consequence of the closeness of symplectic two- form Q. 

The Hamiltonian density, corresponding to Jo, is defined as 

H 1 = 1T u U t + TC v V t - L 

with the result 

#1=2 v2uxx ~ £u Hl = \ ^ 2 ~ ex2 ) Uxx > ( 3 - 9 ) 

where equivalent Hamiltonian densities differ only by a total x-derivative, 
which vanishes in the Hamiltonian W = JJ f^°° H l dxdydz due to appropriate 
boundary conditions at infinity. Indeed, 

1 2 ^ n A 2 ^ ^ 

- X u xx = u + D x \-x u x — xu I •<=>■ u. 

Thus, the mixed heavenly equation in two-component form (12.11) can be 
presented as the Hamiltonian system 

where 5 U and 8 V are Euler-Lagrange operators [6] with respect to u and 
v applied to the Hamiltonian density H\ (they correspond to variational 

derivatives of the Hamiltonian functional J HidV). 

v 
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Bi-Hamiltonian representation of mixed 
heavenly system 



By a theorem of Magri [13J, we can generate second Hamiltonian operator by 
acting with the recursion operator (12. 6p on the Hamiltonian operator ( 13. 8 p 

where Jf 2 , in an explicitly skew-symmetric form, is defined as 

r22 i r /„ „ i i — \ / _ i i 



J i = 9 i - 1 Q- D - — + — D *Q- + v * D v — + — D v v * 



v xx ^XX 



I n 1 1 n \ ( v x u xz u xz v 
- v v D z 1 D,v v - D v 1 D 



'U j xx 'U J xx J \ 'U J xx '^XX ^'XX '^xx 



^x r-, "U>xy Hxy n V x 



+ ( _!L Dz _av + _^ Dz _^ ) (. ; (4.2) 

'^XX '^XX ^^XX 



where we have denoted 



2v 2 1 

Q- = Q = — ( - v x + v z u xy - v y u xz + e) . 

V'xx U>xx 

The proof of Jacobi identity is straightforward and lengthy. The calculations 
can be simplified by using P. Olver's criterion (theorem 7.8 in [6]) formulated 
in terms of functional multivectors. 

We have also performed a straightforward check for compatibility of two 
Hamiltonian operators Jo and J\ using P. Olver's criterion (corollary 7.21 
in his book [6]) and proved that every linear combination a Jo + hJ\ with 
arbitrary constant coefficients a and b satisfies Jacobi identity, i.e. J and J\ 
form Poisson pencil (called Hamiltonian pair in [5]). 

The mixed heavenly flow (12. ip can be generated by Hamiltonian operator 
Ji from the Hamiltonian density 

H = (c - x)vu xx , (4.3) 

where c is a constant, so that the mixed heavenly equation in the two- 
component form (12. ip admits two Hamiltonian representations 

^ -Jo(^) (4.4) 



v t I " u V 5 v Hi J 1 V 5 V H 
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and thus it is a bi-Hamiltonian system. 

We note that we could drop out the term h = cvu xx in the Hamilto- 
nian Hq and set c = in (14. 3p . so that Hq = because the vector of 

variational derivatives of ho belongs to the kernel of J\ 



Ji 



S u h Q 
5 v h 



(4.5) 



According to Magri's theorem, by repeated applications of recursion op- 
erator to the first Hamiltonian operator Jo, we could generate an infinite 
sequence of Hamiltonian operators 



Jn — R n Jo- 



(4.6) 



In particular, for n = 2 we obtain a new Hamiltonian operator J 2 = R 2 Jq 
RJi, which has the following explicitly skew-symmetric form 

DzHD y v z -D,v, 



■h 



^ 2^X i^xyD Z + D z U X y 

UxzDy DyU X z)D ' x 



{{v z D y -v y D z )D- 1 



\ + vtx~ ( V V U xz - V Z U. 

where J| 2 is defined by 



z"yj 



U K^y^xz u z "'xy 



722 

Jn 



(4.7) 



j22 _ V x 



2 - -^T 1 (v y u xz - v z u xy - e)D x + \ v z u xx - -v x u xz j D y 

i^Vy^xx -^z^ ^^D x (Vy1l xz V z U xy £") 



Dy y V Z U XX ^V X U XZ J D z ^ V yU xx ^V x U xy J J ^2 ' 

(4.8) 

Surprisingly enough, the Hamiltonian density, from which the Hamiltonian 
operator J 2 generates the system (12. ip . is proportional to Hi defined by (13. 9|) 



H-i = -eH 1 = 11 — v 2 u xx , 



(4.9) 
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so that the mixed heavenly system admits the three-Hamiltonian represen- 
tation 

The result (14.91) for H_i is derived in a regular way in (15.39P (at c(y, z) = 0) 
in subsection 15.31 (see also (16. 6p ). 

Computing J n in (14.61) for n = 3, 4, . . ., we can generate an infinite series 
of Hamiltonian operators, which shows that the mixed heavenly equation, 
considered in a two-component form, is a multi-Hamiltonian system in the 
above-mentioned sense. 



5 Symmetries and conservation laws for 
mixed heavenly system 

Using the software packages LIEPDE and CRACK by T. Wolf [T7|, run under 
REDUCE 3.8, we have calculated all point and Lie-Backlund second-order 
symmetries of mixed heavenly system (12.11) . For point symmetries, we list 
generators and two-component symmetry characteristics $ = (<p, ip) T , where 
T means transpose: 

Xi = td u + d v , <fi = t, if)! = 1 
X 2 = -d x , ip 2 = u x , if> 2 = v x 
^3 = a (Vi z ) d u, f3a = a(y, z) } ip 3a = (5.1) 
XI = a y {y, z)d z - a z (y, z)d y , ip 4a = a z u y - a y u z , ifj 4 a = a z v y - a y v z 

x 5 = - ^{y9 y + zd g ) + td t + , 

V? 5 = tv - u + ~(yu y + zu z ), if) 5 = tQ + ^(yv y + ZV Z ) 

Xq = td t — xd x — vd v , <£>6 = xu x — tv, ipQ = xv x — tQ — v 
X C ( X)V ) = -c v (x, v)d t + (c - vc v )d u , ip c = c(x, v), ip c = c v (x, v)Q, 

where c(x, v) is an arbitrary smooth solution of the equation 

c xx (x,v) +ec vv (x,v) = (5.2) 
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and we have used the equations of motion (12.11) for eliminating u t and v t . 
Such obvious symmetries as translations in y, z and u, which do not appear 
explicitly in the list (15. ip . can be obtained as simple particular cases of X% 
and X%, while translations in t are obtained from X c at c = —v. 

All second-order Lie-Backlund symmetries [18] modulo point symmetries 
have generators of the form 

X a = a(t, x, v, u x )d u + (a t + a v Q + a Ux v x )d v H , 

(f a = a(t,x, v,u x ), i> a = a t + a v Q + a Ux v x , (5.3) 

where the dots denote an infinite prolongation part and a(t,x,v,u x ) is an 
arbitrary smooth solution of the equations 

O'tx ^■O'vux 0, atv ~\~ a XMx 0, a xx -\- sa vv 0, an -\- &a UxUx 0. 

(5.4) 

Corresponding Lie equations have the form of a second-order flow, due to 
the definition of Q in (12. ip . 

u T = a(t,x,v,u x ), v T = a t + a v Q + a Ux v x (5.5) 

where the "time" r is the group parameter. We note that mixed heavenly 
system itself in (I2.ip is a particular case of (15. 5p at a — v, that obviously 
satisfies conditions (15 .4p . We also note that (<p 6 , ip 6 ) correspond to the partic- 
ular case of second-order symmetries (15. 3p with a = xu x — tv. We point out 
that the symmetry characteristic (ip c , ip c ) in (15. ip together with the condition 
(15. 2p is a particular case of the symmetry (ip a , ip a ) in (15. 3p with a = c(x, v) 
satisfying the conditions (15. 4p . Similarly, we see that symmetry generators 
X\ and X 2 in evolutionary form are also particular cases of the second-order 
symmetry generator X a i t , x ,v,u x ) while X% , X^ v ' z ^ and X 5 are not. 

All second-order Lie-Backlund symmetries can be obtained by taking lin- 
ear combinations of the generators (I5.3P and point symmetries generators 
X 3 " M , X2 {y ' z) and X 5 . 

In table 1 we present commutators of symmetry generators, where the 
commutator [Xj,Xj] is given at the intersection of zth row and jth column. 
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Table 1: Commutators of symmetries of mixed heavenly system. 

y y vf(y> z ) v b (y, z ) y v v v 




















-X, 


x Cv 


Xb v 


x 2 

















-x 2 


-x Cx 


-Xb x 


Y a (y< z ) 








d{a,b) 


2 A 3 

















-y 8(y,z) 











Y a (y< z ) 






d(i,a) 


B(a,b) 


















J£ 9(y,z) 


2 A 4 











x 5 








1 yf 

~2 A 3 


1 Yb 








x c 


—Xb' 


x 6 


x x 


x 2 














x- c 




Xd(x,v) 


~x dv 


Xd x 








-x d 







~ X {b,d) 


X a (t,x,v,u x ) 


—X av 


Xa x 








X a i 


— Xd X( afi ) 


X<£.a,b-» 



We have used here the following shorthand notation: 
<9(a, b) 



d(y,z) 



dyb z (l z hy^ 



CI, c) 0,fC v ®"u, x C-xi 



d(a,b) d(a,b) 

< a, >= — - — - + — -, a = ta t - xa x + u x a Ux - va v 

o(t,v) o{x,u x ) 

a = ya y + za z — 2a, a' = ta t + u x a Ux — a. 



(5.6) 



Let X$ be an evolution form of a symmetry generator with the characteristic 



1> 



Consider the evolution system of PDEs 



v t 



F([u],[v],r) 



(5.7) 



where r = (x, y, z) and [u], [v] denote unknown functions u and v of f together 
with their partial derivatives with respect to the components of r. Let Xp 
be an evolutionary generator of the flow (15. 7p with the characteristic F = 
f 



fj 



. Then X$ generates a symmetry of the flow (5.1) if and only if the 



following commutator relation is satisfied [T9|[20 



X§, Xp 



(5.* 
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where $ 4 denotes partial derivative of $ with respect to t. In particular, 
if $ does not depend explicitly on time t, then these two generators must 



commute: 



Xs>, Xf 



= 0. For mixed heavenly system (12. ip we have F = 

(v, Q) T and then the results of table 1 for Xj, = X v show that equation (15. 8p 
is indeed satisfied for all the generators in the left column of the table, which 
provides an independent test for all of the symmetries of (12. ip . 

5.1 Action of recursion operator on symmetries of the 
mixed heavenly system 

In order to obtain correctly the action of recursion operator R on the space 
of symmetries, we note that the operator D~ x in R should be understood 
as an indefinite integral with respect to x with the "constant" of integration 
C(y,z,t). Indeed, the integrability condition for the recursion relation (ll.4p 
in a two-component notation (12. 5p has the form (p xt = ip x , which implies 
(p t = i> + C(y,z,t) with an arbitrary function C(y,z,t). Therefore, the 
required relation (p t = ip does not follow from (12.51) but may present an 
additional constraint on the t-dependence of C(y,z,t) when R is applied to 

symmetries $ = ( ^ J , while the additive term c(y, z) in C(y, z, t) will still 



be completely arbitrary. In the Lax representation (12.71) and definition (14.11) 
of Ji we have to choose c(y, z) = 0. 

Proceeding to the action of the recursion operator on symmetries, we 
start with the recursion 

<Pi \ _ R ( ¥i\ _ R ( t\_{ D- l u xx \ _ ( u x + C(y, z, t) 



ipi J \ i>\ J V 1 / V Vx ) V Vx 

(5.9) 

but due to the constraint (pt = ip we have Ct = and so C = c(y, z) with an 
arbitrary function c(y,z). Finally we have 

!) = (- + ^*>) = (E) + (£), (5,0, 

where I ^ 3c J = ( Z ^ J . Our next step is to compute 



^3c J V 
^2 J V ^2 J \ V x V ~ £ 



<h\- R ('*) =R («.) = ( ^'(0) ) = ( C(y,z,t) \ 
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using in R the definition of Q from (12. ip . Here the "constant" of integration 
C(y,z,t) could not be neglected since the choice C = will violate our 

constraint (ft = so that ( ^ j will not be a symmetry. Instead, the 

constraint yields C t (y, z, t) = —e so that C = —et + c(y, z) with an arbitrary 
function c(y,z), which implies the result 



(5.12) 



Similarly, we obtain 



At the next step we obtain the first nonlocal symmetry 

\ _ ji ( ¥ia\ _ R ( a z u y - a y u z 



(5.14) 



^4a / \^Aa ) \ 0> z V y - CL y V z 

D~ l {u xx ip Aa + u xz (a y v x + (f Aa) y) - u xy (a z v x + ip Aa ,z)} 

V z {^Aa,y ~ CtyVx) ~ V y ((p Aa ,Z ~ a z V x ) - Q^ A a,x 

with (f4a t y = D y (ip Aa ) and so on, since the expression in curly braces in the 
first row is not a total x- derivative. Another nonlocal symmetry is obtained 
by the action of R on X 5 : 

^ 5 ^-r(^) (5.15) 



V>5 / V ^5 . 

_ 1 ( D- X {u xx {yv y + zv z ) - v x w x + v z w y - v y w z } + etx 
2 V v x (yv y + zv z ) - Qw x + v z w y - v y w z 

where we have denoted w = yu y + zu z — 2u. The action of R on A 6 yields 

A J V ^6 / V U *Q +vv x + ex J \ J 

where we have again used the constraint (p§ t = V>6 on the "constant" of 
integration C(y,z,t). Finally we consider the action of R on the symmetry 
X c in (15.11) with c(x,v) satisfying (15.21) : 

( Pc \ - R ( ¥c \ = R ( c(x,v) \ _ ( D~ 1 {-v x c x + EC V ) 



ip c J \i>c J V c v (x,v)Q J V -c x Q 



(5.17) 
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Equation (15 .2p in the form (c x ) x = —e(c v ) v implies local existence of the 
function f(x,v) such that /„ = c x and f x = —ec v and hence f(x,v) satisfies 
the same equation as c(x, v): f xx +ef vv = 0. The result (j5.17j) . being expressed 
in terms of f(x,v), becomes 

( Pp\-_( D x lD x\f(x,v)] \ _ _ ( f(x,v) \ + f a(y,z)\ ^ ^ 



Vw V SvQ J \ f( x , v )Q J V 

where a(y,z) is the "constant" of integration. If we neglect a(y,z), then 

the result ( 15.181) reads ( % c ) = — ( ^/ | , so that recursion acts on the 

J \ipf J 

solution space of equation (15. 2\i . 

Next we consider action of R on second-order symmetries X a in (15. 3p with 
a(t, x, v, u x ) that satisfies four linear equations (15. 4ft : 

l Pa\ =R f l Pa\ =R { a(t,x,v,u x ) 
4> a J V \cit + a v Q + a Ux v x 



-ea Ux - a x Q + a t v x J ' 

Equations (15. 4p imply the existence of b(t, x, v, u x ) related to a by the equa- 
tions 

b t = -ea Ux , b x = sa v , b v = -a x , b Ux = a t . (5.20) 

As a consequence of ( I5.20p . the potential b(t, x, v, u x ) satisfies the same equa- 
tions (15 .4p as a(t,x,v,u x ). In terms of b, the result ( I5.19P becomes 

4> a J V ht + b vQ + b Ux v x J \b t + b v Q + b Ux v x 

b{t,x,v,u x ) \ f c(y,z) \ ( ifb \. ( <fsc \ /e 91 x 

b t + kq + b UxVx ) + \ o ) \ v& ) v v»3 C ; ' K ] 

so that, up to an arbitrary "constant" of integration c(y, z), the recursion acts 
on the space of second-order symmetries and on solutions of linear system 
fl52D- 

The repeated application of the transformation (I5.20p to a = b takes us 
back to a: b = —ea, so that the action of the recursion on ( I5.2ip results in 
the formula 

R j <£a \ = _ £ ( ^Pa \ _ ( c z u y - c y u z , \ ^ ( d{y, z 




Ipa J ~ \1pa J V C * V V ~ C V Vz ) V 

(5.22) 
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where we have used the relation f 1 5 . 1 3 j) and d(y, z) is an arbitrary "constant" 
of integration. 

The transformed characteristics of Xq, ((p^, ipg) in ( I5.16p . can be obtained 
by the transformation (I5.20p from a(t, x, v, u x ) = ip$ = xu x — tv, which yields 
06 = b — — {u x v + ext). 

The transformation (15.201) . being applied to mixed heavenly system fl2.ll) 
itself, transforms a = v into b = e(x — C), where C is an arbitrary constant, 
so that the transformed flow becomes 

" ^u Ko)=( £( %" C) ) + ( <:< o z) )• (5a) 

where t is the time transformed by the recursion. We note that the action of 
R on (ip c , ip c ) in (15.181) is a particular case of the action of R on second-order 
symmetries in (I5.2ip . 

5.2 Hamiltonian structure of symmetry flows and 
conservation laws 

Hamiltonian operators provide the natural link between commuting sym- 
metries in evolution form [6] and conservation laws (integrals of motion) in 
involution with respect to Poisson brackets. We write Lie equations for sym- 
metries with the two-component characteristics (ipi,ipi) in the Hamiltonian 
form 

where the symmetry group parameter j{ plays the role of ith time for the flow 
( I5.24p and r W — ff /_ °° H l dxdydz is an integral of the motion along the flow 
( 12. ip . with the conserved density H l . The second equality in ( I5.24p is the 
Hamiltonian form of Noether's theorem that gives a relation between sym- 
metries and integrals. We determine conserved densities H l , corresponding 
to known symmetry characteristics ((pi,i/)i), by inverting the relation ( I5.24p 
in the form of inverse Noether theorem 

5 u H l \ = K ( Vi\ = ( u xyViz ~ Uxz<-Piy + ^x^ix + V xx ^i ~ U xx 1pi 

5 v H l J J V 

(5.25) 

where symplectic operator K = Jq 1 is defined in ( 13. 5p . 
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By using (15.251) . we find Hamiltonians for the first four symmetries from 
the list (JED 



Xi : H{ = ( tv - - uj u xx <^=^ H 1 = (tv - xu x ) u xx , 

X 2 : H 2 = vu x u xx , X3 : H^ yz) = a(y, z)vu xx + - (a z u xy - a y u xz ), 

X t : H a(y,z) = ( a zUy ~ d y U z ) ^VU XX + ^(u y U xz - U z U xy ) ^ , (5.26) 

while for the symmetry X5 Hamiltonian does not exist and hence this is 
not a variational symmetry. Equivalent Hamiltonian densities H\ and H 1 
in (15.261) differ by a total x-derivative, which vanishes in the Hamiltonian. 
Hamiltonians for symmetries X c and X 6 are obtained below as specializations 
of Hamiltonians for second-order symmetries. 

Lie equations of second-order symmetries (15. 5p of system (12.11) admit 
Hamiltonian form (15.241) with the Hamiltonian density 

H a = A(t,x,v,u x )u xx -7(t)u, (5.27) 

where 7(t) is an arbitrary function and A is defined in terms of the function 
a(t,x,v,u x ) in (15. 5ft by the relation A v (t, x, v , u x ) = a(t,x,v,u x ). Here A 
should satisfy the equations 

At x — eA VUx = ac(t, x, u x ), At v + A XUx = 0, 
A xx + eA vv =>y(t), A tt + eA UxUx = 5(t, x,u x ) (5.28) 

as a consequence of the corresponding equations (15.41) for the function a with 
restrictions on some "constants" of integration, that follow from (I5.25p . We 
note that at A = v 2 /2 we have a = v and, as a consequence of (I5.28p . 7 = 5 
and the Hamiltonian (15.271) reduces to the Hamiltonian density Hi defined 
in (13.91) . while second-order Lie equations (15.51) reduce to mixed heavenly 
system (12.11) . Therefore, the mixed heavenly system is embedded into the 
hierarchy of second-order flows commuting with it. 

We show now that Hamiltonians for the symmetry flows generated by X c , 
X\, X2 and Xq can be obtained as particular cases of Hamiltonians (15.271) for 
second-order symmetries. The resulting Hamiltonians below are simplified 
by eliminating the terms which are total ^-derivatives. 

Hamiltonian H l of the flow generated by Xi, given in (I5.26p . is a partic- 
ular case of H a in (I5.27P at a = A v = t with A — tv — xu x , that satisfies 
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equations (15.281) with a = 7 = 5 = 0. Hamiltonian H 2 in (I5.26P of the flow 
generated by X 2 is a particular case of H a in ( I5.27P at a = A v = u x with 
A = vu x , that satisfies equations (15. 28|) with a = — e and 7 = 5 = 0. 

Hamiltonian H 6 of the flow generated by Xq in ( 15. ip is obtained from 
the Hamiltonian ( I5.27P of the second-order flow by setting a = A v = ip e = 
xu x — tv which, on account of equations ( I5.28P for A, implies 7 = and 
specifies A and H in ( I5.27P as 

Aq = ^t(ex 2 - v 2 ) + xu x v, H 6 = A 6 (t, x, v, u x )u xx . (5.29) 

Here again we have dropped the terms which are total x-derivatives which 
resulted in the simplification of H 6 and A$ in (j5.29j) . so that the equations 
( I5.28P for A are now satisfied with a = 7 = 5 = 0. 

Hamiltonian of the flow generated by X c in (15. ip is obtained from H a in 
( I5.27P by setting a = A v = c(x, v): 

H c = d(x,v)u xx , where d xx (x, v) + sd vv (x,v) = (5.30) 

and the function d(x, v) is defined in terms of c(x, v) by the relation d v (x, v) = 
c(x, v). Here A = d(x, v) again satisfies equations ( 15. 28ft with a = 7 = 5 = 0. 



5.3 Action of recursion operator on Hamiltonians of 
symmetry flows 

Transformation (I5.10p of the symmetry (<£>i,-0i) implies the following trans- 
formation of Hamiltonian H l 

H 1 = H 2 + H'^ y<z) = vu x u xx + c(y, z)vu xx + - (c z u xy - c y u xz ), (5.31) 

where H^ y , is determined in (15.261) with a(y, z) replaced by c(y, z). 

The action of R on the symmetry (^2^2) m fl5.12j) induces the transfor- 
mation of H 2 

H 2 = -eH 1 + Hl M = -e(tv - xu x )u xx + tf^. (5.32) 

Recursion (I5.13P for the symmetry (<£> 3 , ip s ) i n terms of its Hamiltonian H%, y a 
becomes 

Ha(y,z) — ~Ht(y,z) + ^l{y,z)i (5.33) 
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where H\ yz \ is determined in (I5.26p . For the symmetry (<^ 4 , t/> 4 ), trans- 
formed by R in (I5.14p . the resulting nonlocal symmetry (y5 4 , has a nonlo- 
cal Hamiltonian in the representation (I5.24p with Hamiltonian operator Jo. 
In the next section, we will show that the corresponding nonlocal symme- 
try flow has local Hamiltonian density if 4 with respect to the Hamiltonian 
operator Ji = RJo, as will follow from (16. lip . 

A recursion for Hamiltonian if 6 of the symmetry flow generated by Xq 
will be more convenient to consider at the end of this section as a particular 
case of a recursion for Hamiltonians (I5.27P of second-order symmetries (15 .3p . 

Recursion (I5.18P for the symmetry ((p c ( X)V ), tp c (x,v)) implies the following 
transformation of Hamiltonian if c of this flow defined in (I5.30p 

H c = -H g + Hl {yz) = -g{x,v)u xx + a(y,z)vu xx + -(a z u X y-a y u xz ), (5.34) 

where a{y,z) is an arbitrary function, g(x,v) satisfies the equations g xx + 
eg vv = and g v (x,v) = f(x,v), with f(x,v) being determined in terms of 
c(x,v) by the relations f v = c x , f x = —ec v . 

Recursion formula (I5.2ip for second-order symmetries implies the follow- 
ing recursion for their Hamiltonians: 

H a = H b + Hl M = [B(t, x, v, u x ) + c(y, z)v] u xx 

11 

+ ^ ( c zU xy - c y u xz - 2j(t)), (5.35) 

where B v = b(t,x,v,u x ), b is related to a(t, x, v, u x ) by equations ( 15.201) and 
b satisfies the same equations ( 15. 4ft as a(t,x,v,u x ). 

For the symmetry flow (^gjV'e) i n ( 15. 16ft . obtained by the action of R 
on the flow generated by X$, we take for b = B v in the Hamiltonian ( I5.35P 
the solution b = —{u x v + ext) of equations ( I5.20p with a = xu x — tv. The 
resulting Hamiltonian if 6 for the symmetry flow (15.161) reads 

H 6 = A 6 (t, x, v, u x )u xx + H% M , ^ & = \ i £x<2 ~ v2 ) Ux ~ £xtv > ( 5 - 36 ) 

where equations ( I5.28P for A G = B are satisfied at7 = a = 5 = 0. 

For mixed heavenly system given in Hamiltonian form ( I3.10P with Hamil- 
tonian ifi, defined in f !3. 91) . we apply transformation (I5.20p to a = v to get 
b = B v = e(x — c), where c comes from a constant of integration. Then we 
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determine B in (I5.35p . solving (I5.28P with A n> B, to obtain the transformed 
Hi in the form 

H x = e(x - c)vu xx + Hl (y z) = -eH + H 3 c(y z)} (5.37) 

where Hq is defined in (14. 3ft . The second application of the same transfor- 
mation to (I5.37P takes us back to H\ (modulo the transformed "constant" of 
integration ffL 2 s determined by (I5.33j) ) : 

Hi = \(x 2 - ev 2 )u xx + Hl M = -eH x + H 3 cM (5.38) 
The last equation in (I5.38P is obvious from the alternative formula for H\ in 

(ED. 

For the second Hamiltonian density H = (c — x)vu xx , the transformed 
Hamiltonian coincides with Hi. Hq = H\. This is obvious if one compares 
bi-Hamiltonian representation (14.41) of the mixed heavenly system with the 
second formula in (16. 4p applied to H = H . 

The transformation inverse to (I5.35P : H-i = Hq (with a and b inter- 
changed and b = B v = c — x), yields 

H-i = \(x 2 ~ ™ 2 )u xx + H 3 cM <=► = -e^ + H 3 c(y>z) , (5.39) 

where equivalent Hamiltonian is obtained by dropping a total x-derivative. 

All the flows for the variational symmetries generated by X\, X 2 , X 3 , X 4 , 
-^6, X c ( x ,v)i X a ( t)XjVjUx ) and X 1} X 2 , X 3 , X 6 , X c , X a , where X { = RX i} have 
the Hamiltonian form (15.241) with the local Hamiltonian densities presented 
above, while the nonlocal symmetry generated by X4 has the Hamiltonian 
form with the local Hamiltonian H A 




(5.40) 



with respect to Hamiltonian operator J 1; as will be shown in (16. lip . 

All Hamiltonians of symmetry flows are integrals of mixed heavenly sys- 
tem f[2~T|) . 
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6 Hierarchy and bi-Hamiltonian representa- 
tions for symmetry flows of mixed heavenly 
system 



We know from the work of B. Fuchssteiner and A. S. Fokas (see also 
the survey [20] and references therein) that if a recursion operator has a 
factorized form, as in our case R = JiJq 1 = JiK, and the factors Jo and 
Ji are compatible Hamiltonian operators, then R is a hereditary (Nijenhuis) 
recursion operator. The skew symmetry of Hamiltonian operators j\ = —J\ 
and Jq = — J implies R) = J^Ji = KJi, so that RJ = J\ = JqR) ' . Note 
that in this section, in contrast to sections 15.11 and 15. 3[ we have D~ x = 
dx' (for functions vanishing rapidly at — oo) in the definition (12. 6p of R, 
same as in the definition (14. ip of Ji, so that D~ X D X = 1. 

Now, the Hermitian conjugate hereditary recursion operator 

-D x v x + u xz D y - u xy D z )D x x D X Q - v z D y + v y D z 

*>xxD x 



J{\ — I \ ^x u x T u-xz-LSy — u xy iy z )iy x u x \*t — u z u y ~r v y u z \ 



acting on the vector of variational derivatives of an integral of the flow, again 
yields a vector of variational derivatives of some integral of this flow [19J. 
Therefore, acting with the recursion operator on a variational symmetry flow 

(:)-*(&)■ ^ 

we obtain 



<:)=(s)=*($)=**(£0=*(£* 

(6.3) 

where r and H are the group parameter ( "time" ) and Hamiltonian of the 
transformed symmetry obtained by the action of R on (I6.2p . The transformed 
flow (16. 3 p and its Hamiltonian H have been determined in subsections 15.11 
and I5.3[ respectively, where we now have to skip all arbitrary "constants" of 
integration c(y, z) and a(y, z) for transformed symmetries (<p, ip) T and trans- 
formed Hamiltonians H due to the restricted definition of R given at the 
beginning of this section. Thus, an alternative (to subsection 15. 3j) way of 
transforming a Hamiltonian by the recursion operator is to act by the Her- 
mitian conjugate recursion operator on the vector of variational derivatives 
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of this Hamiltonian 



where both relations follows from (16. 3p . Similarly, 
T f 5 U H\ ( 5 U H \ 
J2 [s v H ) =Jl \5 v H J and SOOn ' 

We now use these remarks for constructing hierarchies of mixed heavenly 
system (12.11) and its symmetry flows, together with bi-Hamiltonian represen- 
tations of the symmetry flows. We start by applying R to this system in 
Hamiltonian form (I3.10P 

(£MS£MS)"(-o°). <«> 

where we have used the result (I5.23p . C is an arbitrary constant, t\ = t is 
the parameter of the group transformed by the recursion and Hi = —eH 
according to (I5.37P with z s skipped, where H is defined in (14. 3p . 

Applying again the recursion operator to the Hamiltonian system (16.51) 
we obtain 

(i)-*(£f:)-*(S:)-*(£S:)-(5)- 

(6.6) 

where ti = t\ and the transformation (I5.38P of Hamiltonian H\ was used: 

H\ = —eH\. This result shows that we are back to the mixed heavenly 
system and so further applications of the recursion operator will not generate 
an infinite hierarchy. We also note that by applying Jo to Ho we will not get 
anything new: 

C:H(SMV)- <-) 

We now apply R to Hamiltonian symmetry flows from (15 .lj) . that com- 
mute with the mixed heavenly flow (I2.ip . to obtain the following results: 

(sMS£MSM&)-(;).m 
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5 U H 2 \ j ( S U H 2 \ _ T ( SuH 1 \ _( 

(6.9) 



vl )- Jl { 5 V H 2 )- Jo { J> ) ~ ~ £Jo V <U? ) -\ I 



U f3 A _ j ( 6 ^ H a(y,z) \ _ j ( 5 U H\ M \ _ ( 5 u H A a 



4a 



(6.10) 



By applying R to the flow of H 4 we obtain the first nonlocal Hamiltonian 
flow in the infinite hierarchy of symmetries of mixed heavenly system: 



where v?4a and ip^ a are defined in (15.141) . The next higher flow of this hierarchy 
can be obtained by applying J 2 to {5 U H 4 ^ 5 V H 4 ) and so on. 

Finally, we consider the action of the recursion operator on the general 
flow of second-order symmetries H 5 . 3 1) with the Hamiltonian (15.27ft 

u t \ _ , ( 6 u H a \ f a(t,x,v,u x ) \ , . 

v t Jo \ 5 v H a J \a t + a v Q + a Ux v x ) " l °' iZj 



Transformation (I5.35P of Hamiltonian (I5.27P under the action of R, H a = H l 
implies the following transformation of the flow 



, , 5 u H a \ j f 5 u H b \ ( bit 

J 0\ S tin ) = J [ , Tjb ) = ( U , \ ^ , A Vgj ) > I 6 ' 13 ) 



5 v H a J u \ 5 V H» J \b t + b v Q + b Ux n 



where b = a is related to a by transformation (I5.20p . with a and b satisfying 

same equations ( 15. 4p . Repeated application of R to (I5.35P yields H a = H b = 
—eH a , because b = —ea, and therefore 

u~ \ = j ( 5 u H a \ = ( 5 u H a \ _ _ e> ( a(t,x,v,u x ) 



v~ J \ s^fia J " \ 5 v H a J \ a t + a v Q + a Ux v a 

(6.14) 

so that we are back to original flow (I6.12p . This is quite similar to the 
behavior of the mixed heavenly system which is a very particular case of this 
general flow. 
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By definition of a hereditary recursion operator, R generates an Abelian 



symmetry algebra out of commuting symmetries. Since 



v a(y,z) v a(y,z) 







this implies 



v a(y,z) ya(y,z) 



0, where X is the symmetry generator ob- 



tained from X by the action of R. Now, X 3 — — ^ 4 



-X. 



a(2/>z) 



for vanish- 



ing constant of integration c(y, z) and therefore 

forward calculation shows that the 
" (04a, 4>4:a) T indeed commute. 



r a(y,z) 
L 4 ) ■ 



ra(y,z) 
L 4 



0. A 



straightforward calculation shows that the flows {u t A,v t i) T = (^Pi a ^4a) T and 



[U t 4 , Vj4 / 



Repeating this reasoning for powers of R applied to the last result, we 
see that the hierarchy of symmetries generated by R from X^ y ' z ^ consists of 
commuting flows. 



7 Two-component form and Lagrangian of 
Husain heavenly equation 

Husain equation (e = +1 in [9]) 

VtyVpz ~ V tz V py + V tt + EVpp = (7.1) 

can be obtained from the mixed heavenly equation ( 11. 3ft by the partial Leg- 
endre transformation in x 

x = -v p , u = v-pv p , p = u x , v(t,p,y,z) =u — xu x . (7.2) 

Up to a change of notation of variables, Eq. (17. ip with e = ±1 is a particular 
case (jl.5p of our general equation (jl.ip admitting partner symmetries. 

To discover its bi-Hamiltonian structure, we have to consider Husain 
equation in a two-component form, without using Ashtekar's Hamiltonian 
formulation of general relativity [2T|[22]. which was a starting point in the 
paper [9] by V. Husain. 

We start with the Lagrangian of Eq. (17. ip in a one-component form 

L = j ( v t + ev l) + 3 v t{ v y v pz ~ v z Vpy), (7.3) 

which for e = +1 is equivalent to the one given in [23] . 
In a two-component evolution form, Eq. (17. ip becomes 

( vt : q (7.4) 

I, It — QzVpy — q y v pz — EVpp . 
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We shall call Eq. (I7.4p Husain system. 

The determining equation and recursion for symmetries of Husain equa- 
tion in a one-component form could easily be obtained by Legendre transfor- 
mation (17.21) from corresponding formulas for mixed heavenly equation and 
then converted to a two-component form, resulting in a Lax representation 
of Husain system (I7.4p . However, the main objects of the Hamiltonian for- 
mulation, like Lagrangian, symplectic two-form, Hamiltonian operators and 
Hamiltonian densities cannot be obtained that way. Therefore, we have to 
undertake an independent study of Husain system along the same lines as we 
have done before for mixed heavenly system. Lagrangian for Husain system 
([73D reads 

L = ^ (2v t q - q 2 + evfj + ^ v t {v y v pz - v z v py ). (7.5) 

Note that the form of the Lagrangian (17.51) is not uniquely defined by one- 
component Lagrangian (17. 3p and it requires some skill in order to arrive at 
the form of Lagrangian that will be suitable for a Hamiltonian form of Husain 
system (17.41) . 



8 Symplectic and Hamiltonian structure of 
Husain system 

Lagrangian (17. 5p yields the canonical momenta 

7T„ = q + - (v y V pz - V z V py ), TT q = 0, (8.1) 

that cannot be inverted for the velocities ut and Vt, and, therefore, Lagrangian 
(17.51) is degenerate. According Dirac's theory of constraints [16], we treat the 
definitions (18. ip as constraints of the second class 

4> v = tt v - q - ^ {v y v pz - v z v py ) =0, <p q = 7T q = 0. (8.2) 

Poisson brackets of the constraints yield the skew-symmetric symplectic ma- 
trix operator 

K = ^ VpyDz ~ VpzDy ~ X Q j . (8.3) 
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The corresponding symplectic two-form Q = J udxdydz of the density 

v 



u — - (v P ydv A dv z — v pz dv A dv y — 2dv A dq) (8.4) 



is closed since the exterior differential of (18.41) is a total divergence (similar 

to d32D) 

duj = — dv z A dv p A dv y , 

which implies vanishing Q at appropriate boundary conditions. Therefore, Q 
and K in (18. 3p is indeed a symplectic form and symplectic operator, respec- 
tively. Hence, its inverse K~ l is a Hamiltonian operator 

Jo = K' 1 = ( 1 V (8.5) 

Indeed, it is explicitly skew-symmetric and Jacobi identity is satisfied as a 
consequence of closeness of the form Q. 

Hamiltonian density, corresponding to Jq, reads 

Hi = ir v v t + n q q t - L = - q 2 - | u^, (8.6) 

so that Husain system (I7.4p is the Hamiltonian system 



where 8 V and 5 q are Euler-Lagrange operators [6] with respect to v and q. 

9 Recursion operator and Lax representation 
for Husain system 

Recursion relation ( II .6p for symmetries of the equation (17. ip . with the change 
of notation u i— >■ v and becomes 

<^p = _ v py<Pz + ft + Uo<P P - (9.1) 
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As before, we introduce two-component symmetry characteristics of Hu- 
sain system determined by the Lie equations (with independent variables not 
transformed by the group) 



Vt ) = ( V 

Qr J V i> 



$, (9.2) 



where r is the group parameter. The symmetry condition is the linear matrix 
equation, the compatibility condition of Lie equations (I9.2p and equations 
flZH) 

T'yZn i(*)=0, (9.3) 

where A is Frechet derivative of the flow (I7.4p 

^ ~~ ( eD* - q z D p D y + q y D p D z , D t + v pz D y - v py D z ) ' 

Then the first row of (19.31) yields <p t = 4>- 

Using the relations ip = ip t , and ip — <Pt for the transformed symmetry, 
we rewrite the recursion relation (19. ip with u>o = and v t = q in the two- 
component form 

¥p = (v pz D y - v py D z )ip + 4>, 

■ip = (q z D y - q y D z - eD p )ip. (9.5) 

In the notation $ = ( ^ j , the recursion relation takes the form $ = -R( ( &) , 
where recursion operator R has the 2x2 matrix form 



R= [ D v\^D y -n py D z ) D-i \ (g6) 
g^D^ - - eD p / 

As in the case of the mixed heavenly system, recursion relation (19. 5 j) for sym- 
metries implies that D" 1 should be understood as an indefinite integral with 
respect to p with the "constant" of integration C(y, z,t) determined by the 
constraint fit = ip up to an arbitrary additive term c(y, z). However, in the 
Lax representation (19. 7p and definition of the second Hamiltonian operator 
Ji = RJ in (HO.ip we should use the restricted definition of D~ x as the 
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definite integral J dp' satisfying the condition D~ X D P = 1 (on assumption 
that all functions vanish at p — — oo). 

The commutator of recursion operator R and operator A of symmetry 
condition f 19 . 3 1) . computed without using the equations of motion, reads 

[R,A}= (9.7) 

D P 1 [( v t ~ q)py D z - ( v t - q)pzDy] , o 

_ (qt ~ QzVpy + q y v pz + ev pp ) y D z - (q t - q z v py + q y v pz + ev pp ) z D y , _ ' 

Thus, on solutions of the system (l7.4j) operators R and A commute and 
therefore R acting on any symmetry $ generates again a symmetry, so that 
R is indeed a recursion operator. Furthermore, vanishing of the commutator 
( 19. 7p reproduces Husain system f)7.4p and hence the operators R and A form 
a Lax pair of the Olver-Ibragimov-Shabat type [TTjlT2"]. 

10 Bi-Hamiltonian representation of Husain 
system 

By a theorem of Magri [T3], second Hamiltonian operator is obtained by 
acting with recursion operator (19. 6p on Hamiltonian operator (18. 5 p 

("0" q ,D t -? y D,-sD,)> < 101 > 

which is explicitly skew-symmetric. The proof of the Jacobi identity has been 
performed by using P. Olver's criterion [6] in terms of functional multivectors. 

We have also made a check for compatibility of the two Hamiltonian 
operators J and J\ by using P. Olver's criterion and proved that every 
linear combination aJo + bJ\ with arbitrary constant coefficients a and h 
satisfies the Jacobi identity, i.e. Jo and J\ form a Poisson pencil (a compatible 
Hamiltonian pair). The flow (I7.4p is generated by J\ from the Hamiltonian 
density 

H = qv p . (10.2) 

Thus, Husain equation in two-component form (I7.4p is a hi- Hamiltonian in- 
tegrable system: 

(:)-*($)-*(&•)■ ^ 
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11 Symmetries and conservation laws for 
Husain system 

Using symmetry package LIEPDE in REDUCE [17] . we have computed all 
generators and two-component characteristics (p, ip of point symmetries of 
Husain system (17.41) : 

X x = d t , (fx = -q, ipx = QyVpz - q z v P y + ev pp 
X 2 = -td t - pd p + qd q , 

V?2 = tq + pv p , 4>2 = q + pq P + t(q z v py - q y v pz - ev pp ) 

X 3 = d p , tp 3 = -v p , ^ 3 = -q p (11.1) 

X 4 = - (yd y + zd z ) + vd v + qd g , 

Vi = v-]^ {yv y + zv z ) , ?p 4 = q-^ (yq y + zq z ) , 
X% = a(y, z)d v , ip ba = a(y, z) } ip 5a = 

Xg = a y (y, z)d z - a z (y, z)d y , (p 6a = a z v y - a y v z , ip 6a = a z q y - a y q z 
X c(t, P ) = (k(t,p)d v + c tt (t,p)d q , y? c = ct(t,p), ip c = c tt (t,p), 

where c(t, p) is an arbitrary smooth solution of the equation 

c tt + ecpp = (11.2) 

and we have used equations of motion (I7.4p for eliminating t> f and qt- Trans- 
lations in y, z and v can be obtained as simple particular cases of X£ and 
Xg, respectively. 

All second-order Lie-Backlund symmetries of Husain system modulo point 
symmetries have the generators 

x f = f(t, p, q, v p )d v + {f t + f q (q z v py - q y v pz - ev pp ) + q p f Vp } d q H , 

<Pf = f(t,p,q,v p ), tyf = f t + fq(q z v py - q y v pz - ev pp ) + q p f Vp , (11.3) 

where f(t,p,q,v p ) satisfies the equations 

ftv p + ef Pq = 0, f pVp - f qt = 0, f qq + ef VpVp = 0, f tt + efpp = 0. (11.4) 

Lie equations corresponding to (111.31) have the form of a second-order flow 

v r = f(t,p, q, Vp), q T = f t + f q (q z v py - q y v pz - ev pp ) + q p f Vp , (11.5) 
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where the "time" r is the group parameter. Husain system is itself a par- 
ticular case of (111.51) at / = q obviously satisfying conditions (lll.4p . so that 
symmetry generator X\ in evolutionary form is a particular case of Xf. Sim- 
ilarly, we note that symmetries X 2 , X 3 and X c u p ) are also particular cases of 
the second-order symmetry -^/(t,p l9 ,«p), while the symmetries X 4 , X^ y ' z ^ and 
Xq v ' z > are not. 

All second-order Lie-Backlund symmetries can be obtained by taking lin- 
ear combinations of generators (111.31) and point symmetry generators X 4 , 
j^-a(y,z) an( j x^ y ' zS> in evolutionary form. In table 2 we present commuta- 



Table 2: Commutators of symmetries of Husain system. 

X\ X2 X 3 X4 X( V ' ^ Xq V ' ^ X c (f )P ) Xg( t ,p,q,v p ) 























Xc t 


Xgt 


x 2 








x 3 











—X c i 


-X~9 


x 3 





-x 3 














x c 


Xg p 


x 4 














1 v f(y,z) 

2 A 5 


1 yHv, z ) 
2 A 6 


-x c 


x s 














d{a,b) 















1 va 
~2 A 5 





J£ S(H,z) 


















d(f,a) 


d{a,b) 















1 va 
~2 A 6 


J£ 9(y,z) 


j£ d(y,z) 








Xd(t, P ) 


~x dt 


Xd> 


~x dp 


x d 











X (d,g) 


Xf(t,p,q,v p ) 


-X ft 




-Xf, 


-Xf 








-X(cJ) 


A '«/,g» 



tors of symmetry generators, where the commutator [X^, Xj] is given at the 
intersection of ith row and jth column. We have used here the following 
shorthand notation: 

d = tc t +pc p -c, g = tg t + pg p - qg q - v p g Vp , a = ya y + za z - 2a 
9 = Q9q + v p g Vp ~ 9, (d, g) = d tt g q + d tp g Vp , 

„, „ d(f,g) d(f,g) . . 
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11.1 Action of recursion operator on symmetries 
of Husain system 

Here we study the action of recursion operator (19. 6p on symmetry character- 
istics of Husain system. We start with the recursion 

<Pi\_ R (<Pi\ =R ( "9 \ =£ (D- l v pp 



tpl J \ 4>l J V % V V* ~ 1z V VV + £V PP J \ % 

l)<V)-il)<l^)> (1L7) 

where the "constant" of integration C(y,z,t) = c(y,z) is time-independent 
due to the constraint <p t = ip. We continue with the formula 

t 2 \ = R ( ^2 \ = R ( tq + pv p 

i>2 J \ 4>2 J \Q + PQp + t(q z v py - q y v pz - ev pp ) 

Dp 1 D p [pq - etv p ] 
p{q z v py - q y v pz - ev pp ) - e(tq p + v p ) 

( Pf=pq~etVp j _j_ | ^ 5c (2/'«) | g\ 

where the first term in the right-hand side of the last equation is a two- 
component characteristic of second-order symmetry (I11.3P with / = pq—etv p , 
that obviously satisfies conditions (111.4jl . Now we compute 



^ ^-R ^ )=r( ~ v ^- I q 



^3 / V ^3 / V ~ < iv J V 1* V PU - % V Pz - eV PP 



/ \ Ipl J \ 4>5c(y,z) 

At the next step, we obtain the first nonlocal symmetry 



11.9) 



4>4 J V^y \ q - | (yq y + zq z ) 

1 f Dp 1 \ v py w z - v pz w y - w t ) 
Q y w z - q z w y + ew p 



11.10) 



32 



where we have denoted w = yv y + zv z — 2v, so that w t = yq y + zq z — 2q. 
Next, we obtain 

&5a(y,z) \ = R ( ^5a(y,z) \ = R ( d(y , z) \ = f d y V z - d z V y 
fp5a(y,z) J V ^5a(y,z) J V / V a V°z ~ a ^V 

C (V, Z ) \ = _ f ^>6a(y,z) \ + f ¥5c(y,z) 

(11.11) 

Another nonlocal symmetry stems out from the transformation of symmetry 

-X-6a(y,z) ■ 

' <ha(y,z) \ = R ( ¥6a(y,z) \ = R ( d z V y - d y V z 

tp6a( y ,z) J V ^ea(y,z) J \ d z q y - d y q z 

D p l [v pz if^ y - VpylfGa^ + ipea} \ (H 12) 

qz<flGa,y — qy(f6a,z ~ £<f6a,p J 

Remaining point symmetry transforms under recursion in the following way: 



\ =R ( <P<%p) \=r( \ = ( D p 1 ^ t ^) 

i>c(t,p) ) V ^<t,P) J \ Ctt(t,p) ) \ -£Ct p {t,p) 

D p l d tp (t,p) \ = f d t (t,p) \ ( Cl (y,z) 
dtt(t,p) J V d u(t,p) J V 

f dm \ ( ^c l{ y,z) \ (1L13) 

where ci(y,z) is a "constant" of integration. Here d(t,p) is related to c(t,p) 
by the transformation d p = Ct, dt = —ec p , so that d(t,p) satisfies the same 
equation (111.2ft as c(t,p): d tt + ed pp = 0. The existence of the potential d(t,p) 
for c(t,p) follows from flll.2|) presented in the form (c t )t = (—£c p ) p . 

Finally, we consider the action of R on second-order symmetries Xf in 
(111.3ft with / = f(t,p, q, v p ) satisfying four linear equations (jl 1.4ft : 

<Pf(t, P ,q,v p ) \ =R f V/(t.M,«p) \=r( ffaP'Q'Vp) 

lpf(t,p,q,v p ) J V Tpffoavp) J \ ft + qtfq + q P fv p 



D p 1 {ft - ev pp f q + q p f Vp ) \ 
+ Qpfq) + Qtfv p J ' 
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where q t = q z v py — q y v pz — ev pp . Equations (111.41) imply the existence of a 
potential g(t,p,q,v p ) for / defined by the relations 

9t = -ef p , g P = ft, g q = f Vp , 9v p = -efq- (11.15) 

As a consequence of (111.151) . g(t,p,q,v p ) satisfies the same equations (111.41) 
as f(t,p, q, v p ). In terms of g, our result fill. 14j) becomes 

( Pj\ =R ( Vf\ = ( D p 1 (9p + v pp9v p + Q P 9q) \ = f D p l D p [g] 

$f J \4>f J V 9t + q P 9v p + q t g q ) V D M 

n M ) + ( C{ Y ] ) = ( Y tm) ) + ( ^ ) > ( n - 16 ) 

D t[9\ J V / V ^9(t,P,q,v P ) J V ^c(y,z) J 

so that, up to an arbitrary "constant" of integration c(y, z), the recursion acts 
on the space of second-order symmetries and on solutions of linear system 
(lll.4p . We note that the second application of this recursion takes us back 
to the original second-order symmetry, up to the factor —e: 

r( P)=s( ? )- ( ^ ) + ( ) (11.17) 

V V>/ J J V ^<y,z) ) V ^d(y,z) ) 

modulo "constants" of integration c(y,z) and d(y,z). 

11.2 Hamiltonian structure of symmetry flows and 
conservation laws for Husain system 

Presenting Lie equations for variational symmetries with the two-component 
characteristics (ipi,ipi) in the Hamiltonian form 

s)-(sMS). (1L18) 

we determine conserved densities H l , corresponding to known symmetry 
characteristics ((pi,ipi), using the inverse Noether theorem 

/ 5 v H l \ = K f <Pi \ = ( V PvViz ~ Vpz^ty ~lj>i\ r u 1q n 
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where the symplectic operator K = Jq" 1 is defined in (18. 3p . Using this relation 
for symmetries from the list flll.ip . we find the corresponding Hamiltonians 
which serve also as conserved densities for Husain system 

X x : H x = - (ev 2 p - g 2 ) = -H u X 2 : H 2 = pv p q +\t(q 2 ~ etj), 

X 3 : H 3 = -qv p = -H , xf y,z) : H^ z) = aq + - (a y v z - a z v y ^v p , 

x < y ,z) , H e = q ^ azVy _ ayVz j _ VyVz (a z v py + a y v pz ), 

X c : H c = c t (t,p)q + ec pp (t,p)v, (11.20) 

where c(t,p) satisfies Eq. flll.2p . For the symmetry X 4 , Hamiltonian does 
not exist and so it is not a variational symmetry. 

Lie equations (111.51) of second-order symmetries of Husain system (17.41) 
can also be presented in the Hamiltonian form (111.181) with the Hamiltonian 
density 

Hf = F(t,p,q,v p ), (11.21) 

where F is defined in terms of the function f(t,p, q, v p ) defined in (111.51) as 
Fq(t,p,q,v p ) = f. On account of equations (111.41) and All. 19j) with H = 
and if = <ff, ip — i>f defined by (111.51) . F can be shown to satisfy the 
equations 

F tVp + eF pq = a(t, p, v p ), F pVp - F qt = 0, 

F qq + eF VpVp = 0, F tt + eF pp = 6(t,p, v p ), (11.22) 

where all the functions in the right-hand sides of Eqs. fll 1 .22j) are arbitrary 
functions of their arguments. We note, in particular, that if H* = F = 
| (q 2 — ev p ) in (lll.2ip . then F satisfies all the conditions ( lll.22p with a = 
0, 5 = and Hamiltonian (111.211) reduces to Hamiltonian (18.61) of Husain 
system, while second-order Lie equations (111.51) with r = t coincide with 
the Husain system (I7.4p . Therefore, Husain system is embedded into the 
hierarchy of second-order flows commuting with it. 

Hamiltonians of the symmetry flows, presented here, serve as conserved 
densities for Husain system. We note that the conservation laws presented in 
this section seem to be different from those given in the paper of V. Husain [9]. 
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11.3 Action of recursion operator on Hamiltonians 
of symmetry flows 

Transformation (I11.7P of the symmetry ((pi,ipx) corresponds to the following 
transformation of Hamiltonian H 1 

H 1 = -eH 3 + Hl {yjZ) = q[ev p + c{y, z)\ + -{c y v z - c z v y )v p , (11.23) 

where c(y, z) is a "constant" of integration in (lll.7p . For Hamiltonian H\ of 
Husain system given in (18. 6p . due to Hi = —H 1 , we have from (11 1.23j) 

Hi = eH 3 - H^ y<z) = -eqv p - H^ z) = -eH - H^ y , z ), (11.24) 

where H is the second Hamiltonian (110.21) in the bi-Hamiltonian represen- 
tation (110.31) of Husain system, that is, Hi = —eH modulo arbitrary "con- 
stant" of integration c(y, z). 

Action of R on the symmetry (y^V^) in (111.81) implies the following 
transformation of the Hamiltonian H 2 

H 2 = + H^ z) = \{q 2 - evl) - etqv p + H%^ (11.25) 

that is, H 2 transforms to the Hamiltonian H? = F in (lll.2ip . where F 
is determined by / = pq — etv p due to the relation F q (t,p,q,v p ) = f and 
equations ( 111.221) for F. It is obvious that the specified / satisfies all the 
equations (111 .4j) . 

Transformation (I11.9P of the symmetry (ip 3 ,ip 3 ) results in the following 
transformation of Hamiltonian H 3 : 

H 3 = -H 1 + Hl M = \{ev 2 v - q 2 ) + Hl M . (11.26) 

Transformation (111. lip of the symmetry (^s,^) implies 

H a(y,z) = ~ H a(y,z) + H c(y,z) = <l( a V V z ~ a ^y) + V y V z (a z V py + d y V pz ) + H\ y z) . 

(11.27) 

Transformation of the symmetry (<y2 6 ,-?/> 6 ) gives rise to nonlocal flow ( 1 1 1 . 1 2 j) 
and hence Hamiltonian H®,s is also nonlocal with respect to the Hamil- 
tonian operator J . In the next section, in Eq. ( 1 1 2 . 6 [) we will see that the 
Hamiltonian of this nonlocal flow will be a local one, notably H%, if con- 
sidered with respect to second Hamiltonian operator Ji. 
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From transformation HI 1 . 13|) of the symmetry (f c (t,p),i , c(t, P )) we deduce 
the transformation of the Hamiltonian H ^'^ 

H c = H d + = q [d t (t,p) + c(y, z)\ + ed pp (t,p)v + -{c y v z - c z v y )v p , 

(11.28) 

where d(t,p) is related to c(t,p) by the equations d p = c t , d t = —ec p and 
satisfies the same equation ( 11 1.2ft as c(t,p): d u + ed pp = 0. 

Finally, transformation ( 111.161) of the general second-order symmetry 
(<ff, ipf), with f(t,p, q, v p ) satisfying equations flll.4j) . results in the following 
transformation of the Hamiltonian H'\ 

= H° + Hl M = G(t,p, q,v p ) + c{y, z)q + -{c y v z - c z v y )v p , (11.29) 

where G q = g and g(t,p,q,v p ) is determined by equations fl 1 1 . 1 5[) for any 
given f(t,p,q,v p ) satisfying ( lll.4p . From (111 .IT)) we note that the repeated 

application of the recursion to H* takes us back to H*\ H* = —eH^ modulo 
"constants" of integration. Since the Hamiltonian H\ of Husain system is 

a particular case of H* with / = q, the same is true for Hi. Hi = —eH\. 
Similarly, the second Hamiltonian H of Husain system can be obtained from 

H* at / = v p , which yields H = Hi and H = —eH . 

12 Hierarchy and bi- Hamiltonian representa- 
tions for symmetry flows of Husain system 

Similar to the beginning of section in this section, in contrast to subsec- 
tions [TO] and [1L31 we have D~ x = dp' (for functions vanishing rapidly 
at — oo) in the definitions (I9.6P and ( 110.1)) of R and Ji, respectively, so that 
D~ X D P = 1. In section El we have noted the identity RJ = Ji = JqB) 
resulting in the relation (16. 4p . which signifies that the action of Ji on varia- 
tional derivatives of the Hamiltonian H can be replaced by the action of Jq 
on variational derivatives of the Hamiltonian H obtained from H by the ac- 
tion of R, in accordance with the formulas derived in subsection 111.31 where 
we now skip all the terms involving arbitrary "constants" of integration. 

We now proceed to use relation ( 16. 4 p for constructing hierarchies of Husain 
system and its symmetry flows, together with bi-Hamiltonian representations 
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of the symmetry flows. Applying R to Husain system in Hamiltonian form 
( 18. 7p . we obtain 



\q tl J \ SqHx J \ 8 q H X J V 5 z H o J \ % J 

(12.1) 

where t% = t is the parameter of the group transformed by R and we have 
used Hi = —eHq due to ( lll.24p modulo H^ y The second application of 
R to Hamiltonian system d 1 2 . 1 [) yields 



v t2 \ h ( 6 Hi \ 6 vHl \ = _ £Jq ( 5 vHl 



q t2 J V SqHi J \ 5„H ) ' V S i H i 




;i2.2) 



due to H\ = —eHx, so that we are back to Husain system and further appli- 
cations of the recursion operator will not generate an infinite hierarchy. We 
also note that applying J to H will not yield anything new 

^ = ( I' ) • (12-3) 



% / V ditto J \ q P 

while Ji applied to Hq is equivalent to Jo applied to Hq = Hi, which again 
results in bi-Hamiltonian representation ( 110.31) of Husain system. 

Next we apply R to Hamiltonian symmetry flows generated by ( lll.ip 
that commute with Husain flow (I7.4p . We skip the flow of H 1 = —Hi, which 
will reproduce ( 112. ip up to a sign, and start with H 2 from ( lll.20p with the 
following results 

v t 2 \ _ t ( <^J^ 2 \ _ t f &vH 



a, J - Jl { 8 q H 2 )- Jo \ 5 q H 2 



pq — etVp 
p{v py q z - v pz q y ) - e(pv p + tq) p 



;i2.4) 



where H 2 = ^ (q 2 — ev 2 ) — etqv p = pH x — etH according to (111.251) . 

Since H 3 = —H , we skip the action of Ji on H 3 . There is no Hamiltonian 
for the symmetry generated by in ( lll.ip . Therefore, we proceed with the 
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action of J\ on H^(y, z) in (111.201) to obtain 

\ - i ( 5 * H kv,*) \ - j ( 5 ^t \__ T ( $vH e a \__( 



QjyV z (X Z Vy 

a y q z - a z g<, 



;i2.s) 



where we have used that H\ = —H® according to (111.271) . 

Applying J\ to the flow of the symmetry (ife a , V'ea), defined in (lll.ip . 
with the Hamiltonian \ given in (lll.20p . we obtain the first nonlocal 
symmetry flow (111.121) in the hierarchy of variational symmetries of Husain 
system 



* )-- h \ 6,Ht )- J °{ SM 



;i2.6) 



because D~ x obviously acts not on a total p-derivative. Since X^ and Xg 
commute, considerations presented at the end of section [H] show that the 
hierarchy of symmetry flows, generated by powers of R from the symmetry 
X£ , consists of commuting flows. 

We skip a discussion of point symmetry X c ^ p ) in (lll.ip . since it is a 
particular case of the second-order symmetry Xf(t >p>q!V _j in (lll.3p . So we 
finally consider the action of R on the Hamiltonian , defined in (lll.2ip 



Vf 

Qt 



Jl (, 5 q Hf )- Jo { 5 q Hf ) ~ Jo V 5 q H* 

g(t,p,q,v p ) 
9t + 9 q {qzVpy - q y v pz - ev pp ) + q p g Vp 



12.7) 



where, according to (lll.29p . W = H 9 with g(t,p,q,v p ) related to / by 

equations (I11.15p . Since fif = —eH*, the second application of R to 
takes us back to the original second-order flow f 1 1 1 . 3 1) with the generating 
function /. We note that the flow (I11.3P is a natural generalization of the 
Husain system. 
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13 Conclusion 



The importance of equations, that admit partner symmetries is that they 
possess nonlocal recursion relations for symmetries of a special form, that 
enable us to obtain noninvariant solutions by symmetry methods. They 
are integrable equations also in a more traditional sense because they admit 
Lax representation together with infinite sets of symmetries and conservation 
laws. Mixed heavenly equation, that combines, up to a point, first and second 
heavenly equations of Plebahski, and Husain heavenly equation are among 
the simplest canonical equations with these properties. We have reformulated 
both equations as two-component evolution systems, which resulted in a 
natural definition of a single matrix recursion operator for each system. This 
operator together with the operator of the symmetry condition forms a Lax 
pair of Olver-Ibragimov-Shabat type for each of these systems. 

By choosing an appropriate Lagrangian, we have discovered symplectic 
and Hamiltonian representations for both mixed heavenly system and Hu- 
sain system. Applying the recursion operator to the Hamiltonian operator, 
we have explicitly generated second Hamiltonian structures for these systems. 
Thus, we have shown that the mixed heavenly equation and Husain equation, 
set in a two-component form, are bi-Hamiltonian systems with compatible 
Hamiltonian structures forming a Poisson pencil. Therefore, they are inte- 
grable Hamiltonian systems also in the sense of Magri. Hamiltonian structure 
relates symmetries and conserved densities which serve as Hamiltonian den- 
sities for symmetry flows. We have determined such Hamiltonian densities 
for all variational point symmetries and generalized second-order symmetries 
for both mixed heavenly system and Husain system. We obtained transfor- 
mations of symmetries and their Hamiltonians under the action of recursion. 

We studied hierarchies of the mixed heavenly system and Husain system. 
A characteristic feature of these systems is that the repeated action of the 
recursion operator on each of these systems takes us back to the original 
system, so that there are only two members in the hierarchy containing each 
of these systems and no infinite hierarchy can be generated from them by 
the recursion. This is a remarkable distinctive feature of the mixed heavenly 
system and Husain system as compared to the second heavenly equation 
of Plebahski and complex Monge-Ampere equation whose bi-Hamiltonian 
structures we studied earlier. However, we discovered an infinite hierarchy 
of Hamiltonian flows commuting with each other and with the considered 
system. Such a hierarchy is generated by one of Lie point symmetries of each 
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system. Higher flows in each hierarchy are nonlocal and we have obtained 
explicitly first nonlocal Hamiltonian flows for both systems. Thus, the mixed 
heavenly system and Husain system possess an infinite number of (mostly 
nonlocal) symmetries and conservation laws which is a customary property of 
integrable systems. We have obtained bi-Hamiltonian representations for the 
flows of all variational point symmetries and all higher second-order symme- 
tries. Further study of higher-order and nonlocal symmetries may provide an 
additional interesting information about the structure of these hierarchies. 

There is also a set of all higher second-order symmetries for each system, 
which have a functional arbitrariness and commute with this system but 
not with each other. The commutators of all the symmetries are presented 
in the tables. This set of second-order symmetries of each system includes 
the original system as a very particular simple case and looks like a natural 
generalization of mixed heavenly and Husain systems. It would be interesting 
to study symmetries and conservation laws of these more general systems of 
second-order Lie equations (15.51) and (111.51) . An important question is if they 
possess partner symmetries, may be under some restrictions. 
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